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Six hundred and seventeeutli Aleettng. 

February 8, 1870. — Monthly Meeting. 

The President in the chair. 

The President called the attention of the Academy to the 
recent decease of Overbeck, at Rome, of the Foreign Honorary 
Members. 

Dr. E. H. Clark made a communication on the results of 
an analysis of one thousand cases of disease in general practice, 
as to the curative action of drugs. 

Dr. Bowditch made a communication, illustrated by a chart, 
on the apparent connection of cloudy days and mortality from 
consumption, for the period from 1811 to 1857, and in this 
vicinity. 

Two papers by Mr. G. W. Hill were presented at the meet- 
ing of December 11, 1869. 

The following problem seems to possess some interest, and I have 
not, in my reading, met with any discussion of it : — 

To determine the elements of the orbit of a planet or satellite, 
which moves in a circle in the plane of the ecliptic, from three obser- 
vations of its direction from the earth, made at equal intervals of time ; 
the positions of the earth and the central body at these times being 
known, but the sum of the masses of the central body and the planet 
or satellite being unknown. ■ 
Or, geometrically stated, — - 

In a plane, given a point as centre and three straight lines, required 
to describe a circle, so that the arcs intercepted between the first and 
second, and the second and third, lines may be equal. 

Let generally B denote the distance of the central body from the 
earth ; 
" " L its longitude as seen from the earth ; 

" " r the radius of the orbit of the planet ; 

" " X its longitude as seen from the earth ; 

" " X i'® longitude as seen from the central body. 

Moreover, employ the subscripts (_i) , („) , Q , to denote the special 
values of the above quantities, which have place respectively at the 
three times of observation in their order. 
VOL. vni. 26 
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If a perpendicular be let fall from the central body on the straight 
line which joins the earth and the body whose orbit is to be deter- 
mined, its length is obviously 

R sin (X — Z) ; 

another expression for the length of the same line is 

r sin (x — ^)' 

Hence for the three times of observation, the three equations 

'•sin (x_i — X_i) =:i?_isin (X.i— Z_i), 
r sin (xo — Xo) = Ba sin (Xo — Zj), 
r sin (xi — Xi) = E^ sin (Xj — Xi). 

But since the orbit is circular, x increases uniformly with the time, 
and consequently xo — X-i = Xi — Xo = "7 suppose. 
Thus the above equations may be written 

»■ sin (xo — 7 — ^-i) = ^-i sin (X_i — Z_i) = a_i, 
rsin(xo — Xo ) = ^o sin (Xq — Z„ ) = Oq. 
r sin (xo + ij— Xi) = B^ sin (Xj — Xj ) = «„ 

which serve to determine the three unknown quantities r, xo, and ij ; 
and it will be noticed that their right-hand members are known quan- 
tities. 

If the sum of the masses of the central body and the body whose 
orbit is sought is denoted by ^, and the common interval of time be- 
tween the observations by t, 

thus, if fi were known, two observations would suffice to determine the 
orbit ; but if /x is not known, ij must be regarded as an independent 
unknown quantity. Hence the necessity for the restriction put at the 
end of the statement of the problem. " Also by this restriction the 
problem is made to depend on the solution of an algebraical equation 
instead of a transcendental one. 

The equations can be simplified by taking two unknown quantities, 
«> and cr, instead of xo and ij, such that 

Xi-f X_i 



= Xo — 



(7 = 1; — 



2 
Xi — X_i 
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and putting 



2 
Then the equations become 

r sin (<» — 0-) = a_i, 

r sin (» -|- 8) = Oq, 

r sin (a -j- 0-) = «!. 
Or, 

r sin a> cos <r = ' — -. 

r sin (w -j- 8) ^ «o» 

r cos <» sin (T = s • 

If r sin 0) and r cos a are eliminated from these equations, and we 
make 

''l "l~ <*— 1 •. n 

-^ : cos a = a = c cos S, 

2ao 



Oil — a. 



- sin 8 = b = c sin /3, 



2ao 

where c may be taken as positive and the quadrant of j3 becomes deter- 
minate, or j3 may be assumed between the limits ± 90°, there will be 
obtained, for the determination of o-, the equation 
sin 2 0- = 2 c sin (<r + ^)' 
The computation of c and ^ may be facilitated by introducing the 
auxiliary quantities h and f, such that 



^ sin f = 



\/\ 



Oq 



k cos f::= -~ — , 
V2ao 
then 

c cos /3 = A cos (45° — f) cos 8. 
c sin ^ = k sin (45° — f) sin 8. 

It is evident that the determination of <r depends on the solution of 
an equation of the fourth degree ; but its value can be very readily 
obtained from the above equation by the tentative process ; and then 
r and by means of the equations 



r sm a 



aok cos (45° — Q 

cos cr ' 
do h sin (45° — () 



r cos a 

sm (T 



and finally xo and >; by means of the relations given above. 
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There is a very simple geometrical construction of the roots of the 
equation in <r. Making cos <t = Xj ^nd sin o- = y, the values of x and 
y are the co-ordinates of the intersections of the curves whose equa- 
tions are 

(X — a) (2^ — b) = a b 

Consequently, if we construct the equilateral hyperbola whose equa- 
tion is 

xy=±h 

and from a point on it, whose co-ordinates are 



x'-- 



y=- 



V^±ab' 

b 
V'TlTb' 



as centre, we describe a circle, whose radius is . , and then draw 

V/±ab 

radii to the points of intersection of the curves, the angles made by these 
radii with the x axis of co-ordinates are the values of er. Since the 
centre of the circle is on the hyperbola, there are at least two intersec- 
tions, and thus the equation in <r has at least two real roots. The geo- 
metrical construction readily affords the condition which a and b must 
satisfy in order that tliere may be four real roots. The condition is, 
that the length of the sti-aight line drawn from the point a, b, on the hy- 
perbola whose equation is 

xy=ab 

normal to the opposite branch, shall be less than unity. The equation 
to the normal which passes through the point x"> y" on this curve, is 

x"(x-x")-!/"(y-f) = o. 

The condition that it passes through the point a, b, gives 

x"(x"-a)-/(/-b) = 0, 
x"/ = ab. 

If we multiply the first of these by x"^, we get 

X"'(x" — a)-ab(ab-bx") = 0, 
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or, rejecting the useless factor x" — 3; 

whence 

X = — Vab^ 
and by interchanging a and b, 

/ =- - ^IFh. 

And thus the length of the normal 

^{x'-^r^-ii/'-^f = [(a + ^■^')' + (b + <^ aTb)']* 

= [at + bt]i 
Consequently, 

if a^ 4" '^^ ^ 1' tli6re will be four real roots ; 

" at -[- b^ = 1, there will be four, and two will be equal ; 

" at -(- by ^ 1, there will be only two real roots. 

We will now show how to arrive at a direct solution of the problem 
by the employment of trigonometric formula. If tan o- is taken for 
the unknown quantity, the equation, on which the solution of the prob- 
lem depends, is 

[c cos /3 tan o- -j- c sin ^f (1 -|- tan* o") = tan* <r, 

or if we put tan o- = X) 

(x + tan^)=(x* + l)=^, 
or, expanded, 

X* + 2 tan ^.X^ + J^ x' + 2 tan ^.x + tan*/3 = 0. 

A quantity n may be assumed, such that this biquadratic shall be 
resolved into the two quadratics 

„ , n sin u cos O 4- /x) I . 

X + 2 — a » X + tan ^ tan /x = 0, 
'^ ' cos ^ cos 2 /i. '^ ' "^ 

„ , „ cos u sin (3 — u) , 

X + 2 o o X + tan ^ cot ;x = 0. 

'^ ' cos ^ cos 2 /* "^ ' '^ 

That this is possible will be evident on multiplying the left-hand 
members of these equations together, for after some reductions easy to 
make, all the coefficients, with the exception of that of xS will be found 
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to be identical with those of the biquadratic ; and consequently ^ is 
determined by the equation 

- Tx I ;. T I o sin 2 a sin O — u) cos (;3 4- u) c' — 1 

tan 13 [tan ^ + cot;.] + 2 cos' ^ cos' 2 / = 3^^^' 

or 

c' sin 2 /3 c' sin 2 ;ii [sin 2 /t — sin 2 ;3] ^ ' 

sin 2 jtt 1 — sin* 2 ^ ' 

or 

sin" 2 /* + (c^ — 1) sin 2 (x — c^ sin 2 /3 = 0. 

That this cubic will always give at least one real value for /», is evi- 
dent on making in the left-hand member sin 2 fj, successively equal to 

— 1, 0, and -[- 1 ; the results obtained are 

— c^ (1 -f- sin 2 p), always negative ; 

— c^ sin 2 /3, negative or positive, according to the sign of sin 2 (3 ; 
-\- c^ {1 — sin 2 p), always positive. 

Moreover, it is plain that there is one real value of /», which makes 
sin 2 n and sin 2 /3 have like signs ; this value we shall adopt. 
Making, according as c^ is greater or less than unity, 

c^ = sec^. y, or c^ = cos^ y ', 

the above cubic is solved by these formulfe (see Chauveuet's Trigo- 
nometry, p. 96), it being necessary to make three diflferent cases. 

Caae I. 

2 sin' V tan y 
tan <b = — , ' '-, 

V27 sin 2 j3 

tan ^ = tan ^, 

2 

sin 2 a = "7= tan y cot 2 ilr. 

Case n. 

2 sin y' tan' y 

sin 4> = — 7=r ' , 

V 27 sin 2 ^ 

tan yjf z= ta,n^ -^, 

2 . , 

sin 2 u = -7— sin y cosec 2 J'- 
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Case III. 

2 sin y tan^ y 

• 2 
sin 2/it = — sin y sin (<^ ± 60°). 

V o 

When (^ is impossible in Case II., the formulae of Case III. must 
be used ; and the upper or lower member of the double sign in the 
second equation must be taken according as sin 2 ^ is positive or nega- 
tive ; in order that sin 2 p. may have the same sign with sin 2 ^. All 
the auxiliary angles <|), •^, and /x may be taken between the limits 
± 90°. Since sin 2 jS sin 2 /* is always positive, tan /S tan \i. and tan 
^ cot ft. are so likewise, since they are respectively equivalent to 

sin 2 sin 2 /[I , sin 2 /3 sin 2 |i 
4 cos^ ^ cos^ ju 4 cos' j3 sin' /i' 

Let us take two auxiliary angles 6 and 8, determined by the equa- 
tions 

tani j3 tani jti cos cos 2 y. 



sin 2 fl ::= ■ 



sin /i cos (;3 -|- ^) 



„ ^ tani j3 coti ju cos (3 cos 2 /x 

cos ju sin O — jn) ' 



or by the equations 



. a ^ — cos 2 
sm 2 ^ = 



fi /sin 2 g 

Y fi) Y sin 2 /i' 



^ cos (/3 + 

. - ., _^ cos 2 u /sin 2 j3 

sm 2 5 = T • .^ s t/ . ^ , 
sm (j3 — j«) y sin 2 jn' 

where the upper or lower of the signs must be taken according as 

^^— ^ in the first and in the second are positive or negative ; 

sm ^ cos ^ r o 

and 2 6 and 2 6 may also be taken within the limits ± 90°. The four 
values of x or tan o- are then 

tan tr =: tani ^ tan^ ^ tan 6, 
tan (T = tani ^ tan* /t cot ^, 
tan a = tani /3 coti ^ tan ^, 
tan tr =i tani /3 coti /t cot 6 . 

If the value of sin 2 5 or of sin 2 does not fall within the limits 
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± 1, it indicates that the two corresponding values of tan o- are imagi- 
nary'. The ambiguity in the determination of o- from its tangent is to 
be removed by taking it in that quadrant which permits the equation 

sin 2 o- = 2 c sin '(o- -j- (3) 
to be satisfied. 

Although all these roots will satisfy the equations with which we 
began this discussion, yet they do not all necessarily belong to the 
problem. The reason of this is, that the three equations are not a 
complete statement of all the conditions of the problem. If we denote 
by A the distance of the body, whose orbit we are determining, from 
the earth, we shall have 

A_i = r cos (xo — 1? — X_i) + i?_i cos {X_i — Z-.i), 
Ao = r cos (xo — Xo) + -^0 cos (Xo — A), 

Ai = r cos (xo + >? — Xi) + -Bi cos (Xi — ii). 

The conditions of the problem demand that A_i, A,, and Aj shall be 
essentially positive. Hence, if any system of values of r, xo and i; ren- 
ders any of these quantities negative, it must be rejected. These re- 
jected solutions really belong to the problem when one or more of the 
quantities X_i, Xo and Xj are increased by 180°. In fact, on referring 
to the equations with which we started, we see they are not altered 
when any one of the quantities X is increased by 180°. The geometri- 
cal statement of the problem is more comprehensive than the applica- 
tion of it to the discovery of the elements of circular orbits. Instead 
of the above criteria for the rejection of solutions not applicable, the 
following, which is simpler, may be used, viz. that x always must He 
in the angle between L -\- 180° and X which is less than 180°. 
This example is added for the sake of illustration : — 
Suppose in the case of Venus revolving about the sun we have these 
data. 

Wash. Mean Time. ^ I. log K 

1869 Jan. 1.0 250° 22' 59".l 281° 24' o4".9 9.9926528. 
« June 15.0 94 37 54. 9 84 33 34. 1 0.0069342 

" Nov. 27.0 292 3 21. 2 245 32 49. 3 9.9939666 

There will be found 

loga_i=9.7048977„ log Op = 9.2497072, log aj = 9.8545925, 
logi =0.5426896, f =324° 41'4."52, 8 = 176°35'15."25, 
log a =9.7678074., log b =9.3111404. 



OP ARTS AND SCIENCES : FEBRUARY 8, 1870. 209 

Constructing the equilateral hyperbola whose equation is xy = — 1, 
and the circle whose radius is 2.89, and the co-ordinates of its centre 
x ■= -\- 1.69, y' = — 0.59, we find the two roots of the equation in o-, 
<r = 7^°, o- = 241^°. In fact, the value of a? + bi = 1.0475 shows 
that the equation has, in this case, but two real roots. Pursuing the 
calculation, 

log c = 9.7928205, ^= 160° 44' 24".60, y' = 51° 38' 20".85. 

Case II. is to be used here. 

</, = — 50° 40' 40".00, r/. = — 37° 56' 3".23, ft = — 34° 30' 27".50 
6 = 14° 49' 46".36, ff is impossible, which confirms the preceding 
statement about the number of real roots ; and the values of <r are 

<r = 7° 23' 36".95 and o- = 241° 37' 18".04. 

If we employ the tentative process with the equation 

sin 2 o- = 2c sin (a- -\- 0), 

we shall get <r = 7° 23' 36".97 and (t — 241° 37' 17".95 ; as these 
values are more accurate, we shall use them. The two solutions are 

a. = 1° 16' 6".99, a> = 197° 31' 54".15, 

log r = 0. 6767422, log r = 9.8624217, 

Xo = 272° 29' 17".14, xo = 108° 45' 4".30, 

t,= 28 13 48 .02, , = 262 27 29 .00. 

On applying the above-mentioned criteria, the first solution is seen 
to be inadmissible, it makes Ao and Aj negative. If both Xq and Xj are 
increased by 180°, the solution will apply. The given example has 
then but one solution. Below we give a comparison between the val- 
ues of the elements of Venus's' orbit as found in this example, and those 
of the " Tables " ; the differences are of course to be attributed to the 
neglect of the eccentricity and inclination of the orbit, and in a smaller 
degree to aberration and perturba'ions. 

From the Example. From the Tables. 

Mean Distance from the sun 0.7284868 0.7233323 

Mean Longitude Jan. 1.0 1869 206° 17' 35".30 204° 57' 28".89 

Mean Motion in Julian Year 2091552".2 2106641".438 
VOL. viir. 27 
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New Method for facilitating the Conversion of Longitudes and Lati- 
tudes of Heavenly Bodies, near the Ecliptic, into Right Ascensions 
and Declinations, and vice versa. 

In the computation of a Lunar Ephemeris, the conversion of the lon- 
gitudes and latitudes into right ascensions and declinations forms no 
inconsiderable part of the work to be done. Prof. Hansen, at the end 
of his " Tables de la Lune," has given some tables, with the view of 
diminishing the amount of labor required in this conversion. 

But their employment seems to me to possess little, if any, advantage 
over the use of the ordinary formulae of spherical trigonometry. I pro- 
pose the following method, which perhaps in a slight degree is more 
ready than that of Prof. Hansen. 

Designating the right ascension, declination, longitude, latitude, and 
the obliquity of the ecliptic respectively by a, S, I, b and «, we have the 
following equations 

sin 8 = cos e sin 5 -|- sin e cos 6 sin I 

. , , sin e . ,, , ,, , sin e . ,, ,, 
= COS e sin 6 -| — - sni {I -\- o) -\ — sm {I — b), 



, I + 90° 
tan^— 



tan 2 


cos — 


f 6-fS 
2 


e — 

cos 


- (A -1- 8) 



2 



The first equation is well known, the second is easily derived from 
the known formula, expressed in the usual notation, 

A B sin (s — c) 
tan — tan — = ^^ -, 

2 2 sm s 

when we remember that, in' considering the triangle, formed by the 
heavenly body and the poles of the equator and ecliptic, A, B, s and 

c are replaced by 90° -f a, 90° — I, 90° + ^ — (^ + ^) and e. 

Suppose we were to tabulate the functions cos f sin A and —^ sin A 

for a certain value of « (as 23° 27' 20" which is nearly its value at 
present), and in small side tables put the variations of these functions 

for increments of 1", 2" 9" in e ; we should have the value of sin 8 

by entering the first table with the argument A = J, and the second 
successively with the arguments A. ■= I -\- b and A = I — b, and 
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adding the results thus obtained, after having corrected them for the 
deviation of the value of € from that adopted in the tables. After 
which the value of S could be obtained from a table of natural sines. 
For the case of the moon, the first function would need tabulation only 
between the limits A =:: 0° and A = 5° 17' ; it might be tabulated for 
every 10". The second would have to be tabulated from 0° to 90°; 
it might be given for every minute of arc. The numbers in these ta- 
bles might be rendered always positive by adding a constant to them ; 
as, for instance, 0.1 to the first function, and 0.2 to the second; and 
thus the addition of the three terms of sin 8 be made easier. 

We should then have to subtract 0.5 from the sum, in order to get 
sin 8 ; or we might prepare a special table, which, with the argument 
0.5 -|- sin 8, should give 8. But by the addition of these constants, the 
extent of the tables would be doubled, as it would be necessary to tab- 
ulate the numbers which correspond to negative values of the ai'gu- 

ments. 

/ _L 90° 
The factor by which tan must be multiplied to obtain tan 

— - — is always positive, and, e being regarded as constant, is a func- 

tion of b -\- B, and, for negative values of S -f- ^) '*^ value is the recip- 
rocal of that which corresponds to positive values of 5 -|- 8. Moreover, 
when J -|- 8 is a tolerably small angle, it does not differ much from 
unity, and varies very uniformly. In the case of the moon b -\- S 
rarely exceeds the limits ± 34°, and the common logarithm of this 
quantity lies between 9.9447979 and 0.0552021 ; and its rate of change 
per minute of arc in J -|- 8 varies only from 262 to 289 units of the 
seventh decimal place. We may, with the better advantage, tabulate 
the function 

log cos — log cos — ~ — , 

for every minute of arc of the argument A from 0° to 34°, with the 

precept that it is to be subtracted from log tan -^I- — when J -)- ^ '* 

a positive angle, but added when 6 -}- 8 is negative. It will be neces- 
sary to append to the table the variation of the function for a change 

I a 

in e. The functions log tan (45° -|- —) and log tan (45° -|- -z) can be 
found from the logarithmic tables, but some labor would be spared had 
we •tables which gave log tan (45° -|- — ) with the argument A both 
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in arc and time ; wliich tables would be useful in many other cases, since 
this function is frequently met with in trigonometric formulae. 

The modifications necessary in applying this method to the inverse 
problem of determining the longitude and latitude from the right as- 
cension and declination are obvious. The variations due to the change 
of the obliquity might perhaps be neglected in using the tables, espe- 
cially in the case of the declination, and computed at the end by means 
of the very simple formulae 

da , . 

:t- =: — tan 3 cos a, 

db 

-y- = sm a. 

. de 

Take this example for illustration : — 

On January 14.0, 1871, G. M. T. we have in the case of the moon, 

I = 206° 40' 35."9 e = 23° 27' 19."81 

b=-\-5 3 16.0 From Tab. I., Arg. 5, +0.0808224 

— 1.7 X (A€=-''>9) 
Z+ 5 = 211 43 51.9 From Tab. II., Arg. /-[- J, —0.1046706 

— 11.7XAf +2 
; — 5 = 201 37 19.9 From Tab. II., Arg. I — b, — 0.0733354 

— 8.2XA6 +2 

8 = — 5 34 37.16 sin 8 —0.0971832 

J-|-8= — 31 21.16 log tan 148° 20' 17."95 9.7900662„ 

a = 13" 46" 19M2 From Tab. III.', to be added, 0.0008223 

+ 0.09 X Ae 



log tan 9" 53" 9'.56 9.7908885„ 

The objection to this method is, that so many arguments l-\-b,l — b, 

J + 8, 45° + — , and a from 45° + — are to be formed ; but this is 

confessedly less fatiguing than the taking of tabular quantities from a 

table. 

It may be allowed to notice here a series, which determines a in 

terms of I, Viz. : — 

, , 2^ € ^ 5 + 8 , 
a = / + - tan - tan — ^ — cos I 

— - tan^ - tan^ "T" sin 2 I 

— |tan«^tan'^4^cos3/ 
8 2 2 
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-f- - tan* - tan' — ^ — sin 4 I 
^- &c 



As tan — tan ■ ~l , in the ease of the moon, is always between the 



2 
31' 



limits ± — , the above series is, for this body, quite convergent. 

€— A 



I add the values of the function log 



cos 



cos 



: + A 



, computed for every 



degree from 0° to 35° of the argument A and for e = 23° 27' 20" 





€ — A 






Change of this 




, ""' 2 






function for an inc. 


A 


log . + A 


A 


A^ 


in e of 1" in units 

of the seventh 

decimal. 




1 


.0000000 
.0015736 


+ 15736 


+ 2 


'+ 0.00 
0.19 


2 


.0031474 


15738 
15744 
15751 
15761 


6 


0.38 


3 


.0047218 


7 


0.57 


4 


.0062969 


10 


0.77 


5 


.00787.30 


12 


0.96 


6 


.0094503 


15773 
15789 


16 


1.15 


7 


.0110292 


17 


1.34 


8 


.0126098 


15806 


20 


1.54 


9 


.0141924 


15826 


23 


1.73 


10 


.0157773 


15849 
15874 


25 


1.92 


11 


.0173647 


28 


2.12 


12 


.0189.549 


15902 


31 


2.31 


13 


.0205482 


15933 


32 


2.50 


14 


.0221447 


15965 


37 


2.70 


15 


.0237449 


16002 


38 


2.89 


16 


.0253489 


16040 


41 


3.09 


17 


.0269570 


16081 


43 


328 


18 


■ .0285694 


16124 
16172 


48 


3.48 


19 


.0301866 


49 


3.68 


20 


.0318087 


16221 
16273 


52 


3.88 


21 


.0334360 


55 


4 08 


22 


.0350688 


16328 


58 


4.28 


23 


.0367074 


16386 
16447 
16511 
16578 


61 


4.48 


24 


.0383.')21 


64 


4.68 


25 


.0400032 


67 


4.88 


S6 


.0416610 


70 


5.08 


27 


.0433258 


16648 


73 


5.29 


28 


.0449979 


16721 
16797 
16877 


76 


5.49 


29 


.0466776 


80 


5.70 


30 


.0483653 


82 


5.90 


31 


.0500612 


I69.i9 


87 


6.11 


32 


.0517658 


17046 


89 


6.32 


33 


.0.'i34793 


17135 
17228 


93 


■ 6.53 


34 


.0552021 


+ 97 


6.74 


35 


0.0569346 


+ 17325 




+ 6.95 



